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1 Abstract 

An improved unified formulation based on the effective field theory is introduced for a spin-1/2 
Ising model with nearest neighbor interactions with arbitrary coordination number z. Present 
, formulation is capable of calculating all the multi-spin correlations systematically in a representa- 
^ I tive manner , as well as its single site counterparts in the system and gives results for the critical 
^ , temperature of the system much better than those of the other works in the literature. The for- 
• mulation can be easily applied to various extensions of s-1/2 Ising models, as long as the system 
contains only the nearest neighbor interactions as spin-spin interactions. Keywords: Ising model 
; Effective field theory ;Correlation functions 
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^ 2 Introduction 
o 

^ 1 Ising model p[] has been one of the most extensively studied models in statistical mechanics and 
' condensed matter physics for a long time. The reason is due to the fact that this simple many-body 
system suits well for investigation of thermal and magnetic properties of various physical systems. 
From the theoretical point of view, this model has been handled by a variety of approximation 
■ techniques and each of them were developed for a demand of obtaining better results, especially for 
determination of critical temperature of the system as accurate as the known exact results. Bethe- 
Peierls approximation (BPA) Bethe lattice approximation (BLA) [3] the effective field theory 
(EFT) ^ the correlated effective field theory (CEFT) ^5, the cluster variation method (CVM) 
1^ ■ the series expansion methods (SE) [7 and expanded Bethe-Peierls approximation (EBPA) [8] are 
Ch ' among the most widely used approximations for the s-1/2 Ising systems. Among these methods, 
O ,. EFT in which the system is treated as a finite cluster is regarded as one of the most powerful 
methods which gives rather accurate results than those of conventional mean field theory which 
considers a many-body system as a single particle system. 

In a typical EFT method ^ , one writes the spin interactions on a finite cluster explicitly by 
, choosing an appropriate Hamiltonian. Contributions coming from the outside of the defined cluster 
' are represented by an effective field, and determination of this effective field plays a key role for 
. solving the problem. Another alternative choice is to work on a larger cluster and to define the 
■"sj" ' interactions only within this larger finite cluster [Tn],[II]. In EFT methods, calculations are often 
. carried on by starting with exact spin identities [12], [13]. These identities include the thermal 
average of some certain functions which depend on the spin variables in the cluster. Hence, in 
addition to the single spin cluster identities there are several works based on two or more spin 
cluster identities in the literature [10], [11]. Whatever form of identities used, in order to perform 
the thermal averages, differential operator ^ and integral operator techniques [Il],[in] are often 
used. However, the accuracy and responsibility of the results obtained by these EFT methods 
costs some complicated mathematical difficulties. Namely, it is widely believed that if one tries 
to treat exactly all the spin-spin correlations emerging on the expansion of the spin identities 
then the problem becomes mathematically intractable. In order to overcome these difficulties, one 
often has to make an approximation on the problem of evaluating the correlation functions which 
means a less accuracy of the results in comparison with the exact ones. Hence, it is no matter 
that whatever the technique is used, the most of the studies in the literature refer to a decoupling 
approximation (DA) which neglects the multi-spin correlations [4| and gives the results identical 
to those of Zernike approximation |16j . However there are small number of works try to handle 
these correlations [5], [9]. Although they have different formulations, these two works gives results 
which are identical to those of BPA for the critical temperatures. 
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Recently, we have introduced an improved formulation, based on differential operator technique 
for investigation of the thermal and magnetic properties of a 3D random field Ising model |17] . 
Our formulation improves the results of the conventional EFT methods by taking into account 
the multi-spin correlations. In this paper, we wish to extend the method given in Ref. |17| for a 
general Ising s-1/2 system with nearest neighbor interactions by proposing a unified formulation 
which systematically includes the multi-spin correlation effects for a given arbitrary coordination 
number z. The general formulation being presented here can be easily adopted to more complicated 
systems such as the models with transverse field. For this purpose, we organized the paper as 
follows: In Sec. |3] we explicitly present the formulations. The numerical results and discussions 
are summarized in Sec. 2] with application of formulation to pure system and bond diluted system 
and finally Sec. [S] contains our conclusions. 



3 General formulation 

The Hamiltonian describing our model is 

= - ^ JijSiSj - '^HiSi (1) 

<i,j> i 

Here Si denotes the z component of the spin variable and it takes the values Si = ±1, = Jji > 
is the ferromagnetic exchange interaction between spins i and j, and Hi is the external longitudinal 
magnetic field on a site i. The first summation in Eq. ([l} is over the nearest-neighbor pairs of 
spins and the other summation is over the all lattice sites. The quantities Jij, and Hi may be 
given with certain distributions or they can have the same values for all pairs of spins/sites, i.e 
Jij — Hi H . 

As our model system, we consider a regular lattice which has N identical spins arranged with 
coordination number z. We define a cluster on the lattice which consists of a central spin labeled 
So, and z perimeter spins {ss, S — 1,2, ... z) being the nearest-neighbors of the central spin. The 
nearest-neighbor spins are in an effective field produced by the outer spins, which can be determined 
by the condition that the thermal and configurational average of the central spin is equal to that 
of its nearest-neighbor spins. 

According to the Callen identity [12] , averages of the spin variables in the cluster are given by 

= tanh (2) 

where /3 = l/ksT, ks is Boltzmann constant and T is the temperature, Ei is the part of the 
Hamiltonian in Eq. ([T]) which includes all contributions associated with the site i and {fi} is an 
arbitrary function which is independent of the site i. The inner average bracket (which has no 
subscript) stands for thermal average and the outer one (which has subscript r) is for configurational 
average which is necessary for including the effect of random bond and random field distributions. 

At this stage, Eq. ([2]) is exact but it is not possible to perform the averages without making 
an approximation due to the large number of degrees of freedom of the system. One coarse 
approximation is well known mean field approximation. In this approximation (for the system 
homogenous distributed bonds, i.e. = J for all i,j and zero magnetic field) the spin Si interacts 
with the field Jz (s) where (s) — (si) , i — 0, 1, . . . , z. Thus, in this approximation magnetization 
is given by 

(s) = tanh(/3Jz(s)) (3) 
On the other hand within the EFT, average of the central spin is obtained from Eq. ([2]) with 



Eq = -^Josss - Hq as 



(5=1 



(({/o}so)), = ( ( {/o} tanh I pTjosss + /3i?o ) ) ) • (4) 



(5=1 




By choosing {/q} = 1, we can get the thermal and configurational average of a central spin from 
Eq. In a similar manner, perimeter spin in the cluster, interacts with central spin and with an 
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effective field wliicli is produced by the (2 — 1) spins outside of the cluster. Hence, the average of the 
perimeter spin can be obtained again from Eq. , but this time with i?^ = ~Jossa~{z — l)h — Hs. 

{{{fs}ss)).r = {{{fs} tanh (/SJq^So + Piz - l)h + ^H,))),^ (5) 

where h is the effective field per spin. By choosing {fs\ — 1, we get the average of the perimeter 
spin from Eq. ([5]) . Since all the sites are equivalent in the model, the value of the effective field h 
for a fixed temperature and certain distributions can be determined from the relation 

(N)). = (6) 

Hence, for a given (/3, h) pair with a certain bond and longitudinal magnetic field distribution, we 
can obtain magnetization as a function of the temperature. Since the effective field is very small in 
the vicinity of critical temperature then the critical temperature can be obtained by letting ft, — > 
in Eq. ([6]) and solving it for /3 — /Sc where /3c = ^/ksTc, Tc is the critical temperature. 

In order to evaluate right hand sides of Eqs. (H]) and ^ we can apply the differential operator 
technique [4] 

exp(aV)/(x) =/(a: + a) (7) 

where V = d/dx is the differential operator and a is any constant. With the help of Eq. ^ we 
get the following equations 

(({/o}so)), = (^^{-^o^n [cosli(JoA-V) + S5sinh(Jo5V)]^^ f{x)U=o (8) 

{{{fs}ss))r = {{{fs} [cosh(Jo5V) + So siuh ( Jq^ V)]))^ fix + {z- l)ft)U=o (9) 



where 



fix) - J dH,PHiH,)gix, H,), gix, H,) = tanh [/? (a; + H,)]. (10) 



Here PniHi) is the probability distribution function of longitudinal magnetic field. 

By expanding the right hand side of Eq. ([5]) with {/o} = 1 we can obtain the thermal and 
configurational average of the central spin. This expansion which is superior to conventional MFA 
takes into account the self-spin correlation identity which is given by 

5^ = 1,5^-1=5. (11) 

for all i where n is positive integer and the expansion produces the multi-spin correlations among 
the perimeter spins for s-1/2 Ising system. In order to avoid dealing with the mathematical 
difficulties originating from these multi-spin correlations, an approximation which is widely used 
in the literature, i.e. DA [ISjcan be used according to 

{{s^S,...Su)), = {{s,)),{{s,))^...{{sk)),. (12) 

However, the multi-spin correlations can be evaluated by using an appropriate {fs\ variable and 
expanding Eq. (|9]). In other words, we can calculate any correlation which includes only perimeter 
spins from Eq. ([9]) by using an appropriate {fs}. Moreover, if we expand the right hand side of 
Eq. (|9]) with a selected {fs} some other correlations appear which are represented in terms of the 
central spin and perimeter spins in the cluster. Similarly, these correlations can also be calculated 
by using Eq. ([8]) and by selecting a suitable {/o}- As a result of this process, we obtain a set of 
linear equations in which all correlations are treated as unknowns. Since this approximation takes 
into account the multi-spin correlations, it gives better results than those obtained by EFT with 
DA. 

Now, in order to obtain a system of linear equations with arbitrary coordination number z, let 
us write Eq. ([5]) as. 



(({/o}.so)). = (({/o}E 



?i=0 



]^S5sinh(Jo5V) J]^ cosh(Jo5V) 

S=l S=n+1 



/(a:)|x=o (13) 
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While writing Eq. ([SJ as Eq. (jl3p we use an approximation that aU n muhi-site correlations 



among the perimeter spins are equal to each other and we represent them with (is 



{{siS2 ■ ■ ■ Sn))r- 1^ Other words we represent 
with one correlation. For example for z = 3,n = 2, the correlations 

{{siS2))r = ((S2S3)),^ = {{siS-i))^ 
o(2) 



possible n-site perimeter spin correlations 



correlations which include central site so i-e. ( ( sos^" 



are identical and represented as = siS2- The same approximation holds also for multi-site 

= {{sosi . . . Sn~i))r represents all 

\ / / r 

n-site correlations which include the central site. 

Now, since all the terms except ss are independent of the thermal average then we can write 
Eq. ([T^ in the form: 



n=0 



(") 
'5 



(14) 



where s^""* = siS2 ■ . ■ s„ and the coefficient An is given by 



An — 



Y[sinh{JosV) Y[ cosh(Jo5V) 

.5=1 5=n+l 



f{x)\x=0 



(15) 



The configurational average can be taken by using any given bond distribution probability function 

Pj {J^J) 



Ar. 



n [dJo^Pj ( Joi)] n Nnh ( Jo^V)] J] [cosh(Jo5V)] 



/(a:)U=o (16) 



5=1 5=1 

The derivation of the multi-site correlations by adopting a suitable {/o} in Eq. (IT4l) requires 
the determination of the identity {^{^{/ols^"''^^ which lies on the right-hand side of Eq. (|14l) . Let 

us call choosing {/q} — Sk in this term as applying Sk to ((sj" // ■ Applying Sk to ((sj" // will 



yield another correlation. Resulting correlation may be ( ( s''" 



,(»+!) 



5 //^ //^ depending 

on the relation between k and n. For determining the resulting correlation we use Eq. (jlip and it 
will be given by 



in) 
SkSg 



,("+!) 
'5 

,("-!) 



k > n 
k < n 



According to Eq. P7)) . if we select {/o} = si in Eq. (ITi)) then we get 

z 

{{sosi)), = A, + {Ao + A2) ((si)),, + Y^An ((4""'^))^ 

n— 3 

Similarly, if we apply S2 to Eq. ((TS)) we obtain 

z 

((S0S1S2)), = (Ai + ^3) {{Sl)), + (Ao +A2+ A4) {{S^S2)), + Y^An ({4"^'^ 



(17) 



(18) 



n— 5 



(19) 



After successive operations, we can get a general expression for the correlation ((soSiS2 • • • Sfc))^ as 
follows 



/ (k) 



fc-1 



ra=0 



n=0 



,("-fc) 
'5 



(20) 



Tl = 2fc + 1 
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Multi-site correlations which appear on the right hand side of Eq. (I^H)) can be obtained from 
Eq. Writing Eq. © for 6=1 and choosing {/i} = S2 in this expression will yield 



{{siS2)),. = B^{{s,)),. + B2{{sos^))^ (21) 

where 

Bi = JdJoiPj{Joi)cosh{JoiV)f{x + j)\x=o 



B2 = J dJoiPj ( Joi) sinh ( JoiV)/(x + 7)^=0 



(221 



Applying S3 to Eq. (1211) will yield the three site correlation which include only perimeter spins in 
the cluster: 

((S1S2S3)), = Bi ((S1S2)),, + B2 ((S0S1S2)), (23) 
By successive derivations we find 

which is the correlation expression for k perimeter spins. One other possible choice for deriving 
correlations which include central site is to apply so in Eq. (IMl) [9i . As seen in Ref. [9J this 
procedure will produce the results of BPA. By taking into account Eq. (|lll) we can write these 
correlations as 

((.«4")),-«.((-r")),+B.((4'-')>, (^'^j 

Thus we can get a system of linear equations by labeling the correlations; 



x^+k+i = {{sasf^)) , fc = 0,1,2,...,2: 



(26) 



By rewriting Eq. ([20]) according to Eq. (|26l) we get 

z-2k 

X:^+k+l = Xk-lTk-l+XkCk+ '^A2k+nXk+n, fc = 0, 1,...,Z (27) 

where 

k k 
Ck — ^^2n, Tk = ^A2„+l. (28) 
n=0 Ti=0 

Similarly for Eq. ([Ml) with Eq. (gSl) we have 

Xk = BiXk-i + B2X^+k, k = l,2,...,z (29) 
and for Eq. ^ with Eq. ^ 

Xz+k+i = BiXz+k + B2Xk-i, k = l,2,...,z (30) 

Note that x-i has to be regarded as in Eq. (P7)) . 

Finally, we construct the system of linear equations based on Eqs. (|27p and ((29)) with the 
coefficients given in Eqs. (IT5|) . (P^ and (l^5|) . By solving this system with the condition Eq. (jB]) (i.e 
xi = Xz+i) we can get all representative correlations defined on the cluster and hence, we obtain 
better results than DA for the critical temperatures of Ising systems with coordination number z. 

Another alternative choice for constructing linear equation system is forming it from Eqs. (j27p 
for only k = 0, ((29|) and (|30l) . The results obtained by forming equation system in this way 
corresponds to results in [9]. These are slightly different from results of DA. Instead of this we 
form our equation system from Eqs. (|27p and (22) • This will give better results than DA and 
results in [S] as seen in Section 2] This may be due to the fact that the equation system which 
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formed from Eqs. (P7|) and ([29]) includes smaller number of terms which have effective field than 
the system formed from Eqs. (H?]) for only fc = 0, (HH) and ((50)) . 

For completeness of the work, let us make discussions on the DA. The equation of state corre- 
sponding to DA can be obtained from Eq. p7| for k = and by letting xj ~ as part of Eq. 
(fT2|) where m stands for ((so))^- 

z 

Ao + ^A„to" - m = (31) 

Solving (|3T|) with coefficients given in Eq. ([T6l) will give the magnetization within the DA. The 
critical temperatures can be obtained by letting m —> in (|31]l and solving it for f3 — Bn\18\. 



4 Results and Discussion 

Since the aim of this paper is to introduce an improved formulation within the EFT which is 
capable of calculating the multi-site correlations for any coordination number z, we give only few 
results in comparison with DA and EFT results. Hence, we wish to present the results of the 
formulation presented here only for a pure system with zero magnetic filed and for a bond diluted 
system (with bearing in mind that the full discussion of this problem is out of the scope of this 
paper) in comparison with DA and EFT results. The results of random field Ising model on three 
dimensional lattices can be found in our earlier work |17) . 



4.1 Pure System with Zero Magnetic Field 

This is the simplest s-1/2 Ising system. The distribution functions are delta functions as Pj (Jij) = 
S {Jij — J) for all i,j pairs of spins and Ph {Hi) ~ 5 {Hi — H) where H — ioi all i. Under these 
distributions, the coefficients of our linear equation system which are given in Eqs. (fTC|) and (|22p 
gets the form 

An = f Mcosh^-"(JV)sinh"(JV)/(a;)U=o 

Bi = cosh(JV)/(x + 7)U^o ^ ' 

B2 - sinh(JV)/(x + 7)| 

x=0 

In Eq. ([5^ . if we write the hyperbolic functions cosh ( JV) and sinh ( JV) in terms of exponential 
expressions exp ( JV) and exp (— JV), and use the binomial expansion with the differential operator 
technique in Eq. ^ then we can obtain more suitable forms of the coefficients An , Bi , -62 as 

^. - *(;)£i:(';")('.')™-^^-^'>''i 

\ / r=Os=0 ^ / \ / (33) 

Bi = I [/(J + 7) + /(-J + 7)] 
B2 = l[f{J + l)~f{-J + l)] 

By solving the system of linear equations defined by Eqs. (P7|) and for an arbitrary 
coordination number z with the coefficients given in Eqs. ([^5]) and ([55]) we get the temperature 
dependence of the whole representative correlations defined in the system, as well as the critical 
temperature as a function of z by taking the limit 7 in the coefficients. 

In Table ITTl we compare the critical temperature kBTc/J of the system obtained within the 
present work with those of the other methods in the literature for various coordination numbers. 
As seen in Table 14.11 our formulation gives the best approximated values for the numerical values 
of the critical temperatures when compared with MFA, EFT 9 and DA [18], which originates from 
the consideration of the multi-spin correlations in the system. The exact results listed in Table FOI 
are the results of high temperature series expansion method [19 which has been regarded as the 
exact results in the literature except the case z = 4 which has analytical result [20] . 

The improvement of this formulation is not limited with the better critical temperatures when 
compared with other approximations. The computability of correlations allows us to determine 
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some of the thermodynamie funetions of the system, such as the internal energy and the specific 
heat. 

In Fig. [l] we can see the temperature dependence of the specific heat curves which are obtained 
from EFT and our formulation in comparison with Monte Carlo (MC) simulation for the square 
lattice {z = 4). In MC simulation 100 x 100 lattice and standard Metropolis Algorithm was used. 
Since MC simulation gives behavior for specific heat close to real one, we can conclude from Fig. 
[T]that our formalism gives more accurate behavior of the specific heat than EFT. This means that 
formulation presented here handles two site correlations more accurate than EFT and this results 
in a more accurate behavior of internal energy with temperature, as well as the behavior of specific 
heat with temperature. 
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Figure 1: Variation of specific heat of square lattice with temperature for pure system under zero 
magnetic field. While the dotted curve shows the MC results, the red curve is EFT result and the 
blue one is our formulation's result. 

The temperature dependencies of magnetization, internal energy and specific heat for a simple 
cubic lattice {z — 6), are shown in Fig. [3] , in comparison with DA (4j and EFT [S] as a limiting 
case c = 1 of bond diluted system. 



Table 1: The critical temperatures of pure system with zero magnetic 
field obtained by several approximations as well as exact results and 
the results of the present work. 



Lattice 


MFA 


DA [18] 


EFT [9] 


Present Work 


Exact [T9][20] 


3 


3.0 


2.104 


1.821 


1.504 


1.519 


4 


4.0 


3.090 


2.885 


2.536 


2.269 


6 


6.0 


5.073 


4.933 


4.527 


4.511 


8 


8.0 


7.061 


6.952 


6.516 


6.353 


12 


12.0 


11.045 


10.970 


10.499 


9.795 



4.2 Bond Diluted System 

Let us treat bond diluted system with zero magnetic field. The bond distribution function is given 

by 

Pj {Jrj) = cS (Jy - J) + (1 - c) (5 ( J,j) (34) 
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and it distribute bonds randomly between lattice sites to be c percentage of bonds are closed and 
remaining 1 — c percentage of bonds are open i.e. c is the concentration of closed bonds in the 
lattice. 

Calculating the coefficients given in Eqs. ([16]) and p2|) by using Eq. (p4l) yields 

An = ^ ^ [ccosh(JV) -c+l]"""c"sinh"(JV)/(a;)U=o 

Si = [ccosh(JV)-c+l]/(a; + 7)U^o ^^^^ 
B2 = csinh(JV)/(a; + 7)U=o 

By applying the procedure used between Eqs. ([5^ and (155]) to Eq. (|^ gives 

Bi = f [/(J + 7) + ./(-^ + 7)] + (l-c)/(7) 
B2 = f [/(J + 7)-/(-^ + 7)] 

(36) 

Solving Eqs. (I?fl) and (^5)) with coefficients given in Eqs. and (^5)) gives the correlations 
and magnetization as a function of temperature. Again, the critical temperature of the system for 
a given set of Hamiltonian parameters can be obtained by letting 7—^0. In this way the phase 
diagrams and the representative correlations of bond diluted s-1/2 Ising system with arbitrary 
coordination number z can be obtained. 

The phase diagrams in (kBTc/J — c) plane for three dimensional lattices with coordination 
numbers z = 6 (simple cubic lattice), 2 = 8 (body centered cubic lattice) and z = 12 (face centered 
cubic lattice) can be seen in Fig. [2] As seen in this figure, as the bond concentration decreases 
then the critical temperature values decrease gradually and fall to zero, as expected. As seen in 
Fig. [5] our formulation gives lower values for critical temperatures at all concentration values with 
respect to other two approximations. It is well known fact that for concentration values which are 
below a certain c*, the system exhibits no ordered phase at all. This specific c* value is called 
the critical bond concentration value. Since the formulation presented here gives lower critical 
temperature values at all concentration values with respect to other EFT formulations, it will give 
also higher critical bond concentration values. The critical bond concentration values of different 
lattices can be seen in Table l4?2| in comparison with the other two approximations. 

The internal energy and specific heat curves as function of Hamiltonian parameters and tem- 
perature can be easily obtained within the present formulation, as well as better results for critical 
temperatures since the formulation is capable of calculating the multi-site correlations as well as its 
single site counterparts in a representative manner. This can be seen in Fig. [3] in comparison with 
DA and EFT. As discussed above, since the DA neglects all multi-site correlations it will give zero 
internal energy just above the critical temperature of the system which is physically impossible. 
On the other hand EFT gives more reasonable results than DA. 

In Fig. [3l other than critical temperature values, clear distinction stands out between behaviors 
of internal energy just above the critical temperature. EFT gives energy values more close to zero 
for the temperatures T > than our formulation. Also the difference of the behaviors of energy 
just above the critical temperature gives rise to a difference between specific heat behaviors after 
Tc between EFT and formulation presented here. 

The other distinction shows itself in the ground state values of magnetization and internal 
energy when c values get closer to the critical concentration value c*. For these concentration 
values, our formalism gives lower values than EFT and DA for the ground state magnetization and 
internal energy. This result is expected since, while the c starting to decrease from c = 1, the low 
temperature value of magnetization gradually decreases until c reaches critical bond concentration 
value and our formalism gives higher critical values for bond concentration . The variation of the 
ground state values of this thermodynamic functions can be seen in Fig. |4]for simple cubic lattice 
(^ = 6). ^ ^ 

As seen in Fig. |4l variation of the ground state values of magnetization and internal energy 
with bond concentration c for simple cubic lattice shows no significant difference for high c values, 
i.e. where bond diluted system is close to pure system. In contrast to this, the difference between 
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our formulation and the results of EFT and DA for for ground state values of magnetization and 
internal energy shows itself for lower c values. 

Table 2: The critical bond concentration val- 
ues of bond diluted system with zero magnetic 
field obtained by DA and EFT, as well as the 
results of the present work. 



Lattice 


DA[TH] 


EFT [5] 


Present Work 


3 


0.5575 


0.6623 


0.7622 


4 


0.4284 


0.4774 


0.5470 


6 


0.2929 


0.3095 


0.3458 


8 


0.2224 


0.2303 


0.2520 


12 


0.1504 


0.1528 


0.1633 




c c c 



Figure 2: Variation of critical temperature for simple cubic lattice (z = 6), body centered cubic 
lattice (z = 8) and face centered cubic lattice (z = 12) with bond concentration. 



5 Conclusion 

In this work, we present a general formulation of s-1/2 nearest neighbor Ising system with ar- 
bitrary coordination number z. The superiority of this formulation lies under its capability of 
calculating correlations in a representative manner and this advantage shows itself in the results of 
critical temperatures and the variation of thermodynamic functions with temperature such as the 
magnetization and the internal energy. Formulation covers some quenched disorder effects since 
derivation starts with Hamiltonian ([T]) which include bond disorder. 

Disorder effects are important in material science since disorder (like bond dilution) induce 
important macroscopic effects in materials. Thus it is important to obtain critical values (e.g. 
critical temperatures, critical bond concentrations ) as well as variation of order parameter or 
some other thermodynamic functions (e.g. specific heat) with temperature as much as possible to 
the exact ones. It is a well known fact that it is impossible to obtain exact results for systems 
with disorder in most cases. On the other hand MC or similar simulation algorithms give accurate 
results for these systems but with some computational cost. 
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Figure 3: Variation of magnetization (m), internal energy (U) and specific heat (C) per spin of 
simple cubic lattice with temperature for some selected bond concentration values. The continuous 
lines are the results of the presented formulation, dashed lines are the result of EFT and the dotted 
lines are of DA. 



On the other hand, as mentioned in |21| , |22| some diluted antiferromagnets in uniform external 
magnetic field corresponds to a ferromagnet in a random external magnetic field. Then in some 
cases, one can obtain the behavior of more complex systems by solving s-1/2 Ising model or it's 
variants like random field distributed system. 

This work is not the first attempt to handle the correlations which appears when expanding 
exact spin identities like Although the most of works deal with critical behavior of spin systems 
within the framework of EFT, these works are based on DA which means that neglecting all multi- 
site correlations, there are some works handling these correlations [5, [9]. These works give results 
for critical temperatures as BPA. The importance of calculating these correlations are two fold. 
Firstly one can obtain more accurate critical values about the system and secondly one can obtain 
reasonable values for thermodynamic functions which are obtained from these correlations such as 
internal energy and the specific heat. 

However, although this work is not the first attempt to handle correlations, the formulation 
presented here treat these correlations in a different way which give rise to more accurate results 
than small number of works related to it. Beyond that, we believe that it is important to obtain 
general formulation which covers arbitrary lattice and arbitrary Hamiltonian as long as it includes 
nearest neighbor interaction as a spin-spin interaction. We hope that the formulation and results 
obtained in this work may be beneficial form both theoretical and experimental point of view. 



10 




0.9 0.8 0.7 0.6 0.5 0.4 0.3 



-0.2 



-2.0 







■ DA 




'. EFT J 




. Present ) 




Work ff 
/ 






z=6 



1.0 



0.8 



0.6 0.4 

c 



0.2 



0.0 



Figure 4: Variation of ground state magnetization and internal energy of simple cubic lattice 
(z = 6) with bond concentration. The calculations have been done at the temperature T = 0.05. 
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